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Abstract 

We consider the scattering of the brane fields due to t-channel mas- 
sive graviton exchanges in the Randah-Sundrum model. The eikonal 
amplitude is analytically calculated and both differential and total 
neutrino- nucleon cross sections are estimated. The event rate of quasi- 
horizontal air showers induced by cosmic neutrinos, which can be de- 
tected at the Pierre Auger Observatory, is presented for two different 
fluxes of cosmogenic neutrinos. 



1 Introduction 

The detection of air showers induced by uhra-high energy neutrinos may help 
to solve many important problems, such as propagation of cosmic neutrinos 
to the Earth and their interactions with the nucleons at energies around tens 
(hundreds) of TeV. In this energy region, neutrino-nucleon interactions may 
be strong due to a new physics. There is a large class of models in a space- 
time with extra spacial dimensions which result in a new TeV phenomenology. 
In the present paper we will consider an approach with non-factorizable met- 
ric proposed in Refs. [H Ej, and study a scattering of the SM fields in this 
scenario. 
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The RS model ^1 12] is a model of gravity in a slice of a 5-dimensional 
Anti-de-Sitter space (AdSs) with a single extra dimension compactified to 
the orbifold /Z2. The metric is of the form: 

ds^ = ^^1^ ^2;^ + dy^. (1) 

Here y = r^O (0 ^ 6' ^ vr), Tc is a "radius" of extra dimension, and parameter 
K defines the scalar (negative) curvature of the space. 

From a 4-dimensional action one can derive the relation: 

Ml, = — (1 - e-^-) ^ ^, (2) 

which means that k ~ Mpi ~ M, with M being a 5-dimensional Planck 
scale. 

We will consider so-called RSI model PP which has two 3-dimensional 
branes with equal but opposite in sign tensions which are located at the 
point y = TTTc (called the TeV brane) and at the point y = (referred to as 
the Planck brane). All SM fields are constrained to the TeV brane, while the 
gravity propagates in the bulk (all spacial dimensions). 

From the point of view of an observer located on the TeV brane, there 
exists an infinite number of graviton Kaluza-Klein (KK) excitations with 
masses: 

m„ = x„Ke"'"''"", n=l,2..., (3) 
where Xn are zeros of the Bessel function Ji(x):^ 

Ji{xn)=0, n = l,2... (4) 

By using a linear expansion of the metric, one can derive the interaction 
Lagrangian: 



1 1 °^ 
r- —T 7,(0) Lr \^/)(") f^^l 



Mpi A, 



where 



A^ = Mpie--^''^ (6) 

is a physical scale on the TeV brane. It can be chosen as small as 1 TeV for 
a thick slice of the AdSs, Vc — 12/ n ^ 60 Ipi. We see from (jSj) that couplings 

iThe first four values of a;„ are 3.83, 7.02, 10.17, and 13.32. 
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of all massive states are only suppressed by A~^, while the zero mode couples 
with usual strength defined by the reduced Planck mass Mpi = Mpi/^/Sn. 

The main phenomenological parameters of the model are the scale 
and the ratio 

The present experimental data together with theoretical bounds on the cur- 
vature of the AdSs restrict an allowed region for the variable /i (see, for 
instance, Fig. [T] taken from Ref. 

0.01 </i< 0.1. (8) 

The allowed value of fi is restricted by so-called naturalness and requiring a 5- 
dimensional curvature to be small enough to consider a linearized gravity on 
the brane. Thus, the lightest masses of the KK graviton modes, m„ = Xnfi-^-K, 
are of the order of 1 TeV. 

Our paper has the following structure. In the next Section we consider 
interactions of the SM fields on the brane in the RS model induced by an 
exchange of massive gravitons. The eikonal amplitude is calculated and an 
elastic cross section for different sets of the RS parameters and/or invari- 
ant energy is estimated. In Sec. 3 we use these results to study the scat- 
tering of ultra-high energy cosmic neutrino off the atmospheric nucleons. 
The neutrino-nucleon cross section is calculated and an event rate of quasi- 
horizontal neutrino events expected at the Auger Observatory is presented. 
Our conclusions and discussions is a topic of Sect. 4. 



2 Eikonal amplitude in the RS model 

In what follows, we will employ the zero width approximation for the graviton 
KK resonances. The Born amplitude corresponding to t-channel exchange 
looks like (both the massless mode and KK gravitons contribute): 



oo 



= — + ^E3^. (9) 

The sum in Q converges very rapidly in n, since x„ = vr (n+l/4)+0(n^^) [?]. 
We consider the scattering of two particles living on the TeV brane. Thus, 
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in Eq. (jU)) t means a 4-dimensional momentum transfer which is well-defined 
and conserved. 

Let us underline that in Q we sum spin-two particles with difi^erent KK 
numbers n (non-reggeized KK gravitons). In more general approach, one 
should sum Regge trajectories an{t) which are numerated by n (KK-charged 
gravireggeons) . For the ADD model, it was done in Refs. jT3]. The results 
of Refs. [ini [131 '^^^ t)e reproduced in the limit ari(^) ^2. As for the RS 
model, results on a gravireggeon contribution to the eikonal amplitude will 
be presented in a forthcoming paper ^14j . 

Generally, the massive KK states may decay to a pair of SM particles. 
The partial widths are proportional to m^/A^, where m„ is the mass of the 
resonance. In particular, the partial decay widths to massless gauge bosons, 
fermions, and a pair of Higgs are^ 



ml 



ml 



1607rA2' 

3 

r(/,W ^ HH) = (10) 
^ ^ 4807rA2 ^ ^ 

Here Ny = 1(8) for photons and electroweak bosons (gluons), Nf = 1(3) for 
a lepton (quark pair) mode, and a = 1/2 for identical particles. Then for 
the total width of the massive KK graviton in the RS model, r„, we get the 
estimate (see also 0): 

~m„(0.31/ix„)2. (11) 

Since the sum which we are interested in converges very rapidly in n (see 
a comment after Eq. 0), we conclude from (fTTj) and (jSj) that effectively 
r„/2 < m„. 

The sum ® can be calculated analytically by the use of the formula [7] : 



oo ^ 



^These expressions can be obtained by a replacement Mpf A^^ in corresponding 
formulae derived for large extra dimensions in Ref. [HJ- We have also neglected masses of 
the SM particles, since TOg^ <C m„. 
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where z„ 
obtain: 



[n = 1, 2 . . .) are zeros of the function z ^ Jy{z). As a result, we 

^2 1 hiv) 



AB{s,t) = i^— + 



-t ' 2fiAl^hiv) ^^^^ 

Here Ii{z) {i = 1,2) are modified Bessel functions and v = \^^/fiAj^. Taking 
into account properties of Ii{x), we conclude from (fT^ that a contribution 
of the massive graviton modes dominates at large \t\: 

" ' (14) 



AB{s,t) 



|t|>AtA7i 



Note, we would get another asymptotics in t, namely, AB{s,t) ~ if we 
sum only finite number of the massive gravitons. 

As it was shown in Ref. ^U], it is ladder diagrams that makes a leading 
contribution of the KK gravitons to an amplitude and results in the eikonal 
representation for the amphtude (g^ = —t): 

A^'^{s,t) = 2is j d?he''i^ [l - e^'^^"''')] , (15) 

with the eikonal given by 

x(s,&) = ^ / ^^?^o(g&)^B(s,-g'). (16) 

The proper accounting for the massless mode has been presented in 
Ref. The result is the following: 



A''\s,t)=e 



i<p4 



SttGns r(l — iGNs) 



-t T{l + iGNs) 



JXmass{s,b) 



■ (17) 



Here Xmass{s, h) denotes a contribution the massive modes to the eikonal, Gj^ 
is the Newton constant, and 04 is a 4-dimensional (infinite) phase. The first 
term in the RHS of Eq. (|T7j) is well-known 4-dimensional result derived by 
different methods in Refs. ^2] • It is negligible at any conceivable energy and 
momentum transfer, and we can write (up to a phase factor): 



(18) 
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It follows from (|13|) that the eikonal depends on two dimensionless vari- 
ables, s/A^, and 

u = bjjA^, (19) 

and it looks like 



X,mass\S,U) = X,mass\^) 



oo 

/xA^ Stt A^J /i(t;) 



The eikonal ()20|) is very well approximated by the following expression (see 
Appendix for details): 

~ ' ^ ^ exp(-2V3u). (21) 



IGtt mA^ 

At a/s ^ 5Ajr, the eikonal is exponentially small outside the region 

^^""W^vfk'-^- '''' 

At 6 — i> 0, it is proportional to . Thus, we can roughly estimate the 
high-energy behavior of elastic cross section: 

where m\ is a mass of a lightest KK graviton. 

Let us underline that the Froissart-Martin like formula describes the 
contribution of the massive graviton modes. The presence of the massless 
graviton in the theory should result in infinite elastic and total cross sec- 
tions ^H]- However, its contribution can be safely neglected in our further 
calculations. 

We can rewrite Eq. ()18p in the form 



A^'^s, t) ~ Am , f \ duuJa (u^^^^ \l - e^*™-(^'")l . (24) 



Correspondingly, the differential cross section in dimensionless variable 



y = ^ (25) 



6 



is defined by 



dy 



1 



(26) 



and we get the estimate 



daei 



dy y- 



■0 36 (/iA^) 



vr s 




(27) 



It follows from (j^ . that da^i/dy depends only on variable y (|^. 
parameter /i, and the ratio ^/K-j^, in addition to the dimensional factor 
(/iA^r)"^ which defines a magnitude of the cross section. In particular, we 
have daei/dy\y=o = s (/xA^)"^ /(s/A^), and a^i = A^)'"^ g{s/ Al, /i), were 
f{x) and g{x,y) are dimensionless functions defined via the eikonal. 

The results of our calculations with the use of formulae (j2l|l (j26|l . and 
fl2()|l are presented in Fig. [21E1 The curves in Fig. |21 which show an energy 
dependence of the cross section were obtained for A^ = 2 TeV, and fi = 0.05. 
The dependence of the differential cross section on the parameter A^^ at 
s = 2 ■ 10^^ GeV^, fi = 0.1 is presented in Fig. El Next Fig. ^demonstrates 
the dependence of da /dy on the parameter /i at s = 2-10^° GeV^, A,r = 1 TeV. 
Finally, in Fig. the reduced differential cross section (namely, multiplied by 
the factor (/iA^)^) is shown for several sets (A^,/i). 

3 Neutrino-nucleon cross section and neutrino 
induced air showers 

Let us now estimate a neutrino-nucleon differential cross section as a func- 
tion of variable y. The neutrino scatters off quarks and gluons which are 
distributed inside the nucleon. Thus, neutrino-nucleon cross section is pre- 
sented by 



where /j(x, M^) is a distribution of parton i in momentum fraction x, and s = 
xs is an invariant energy of a partonic subprocess. The partonic differential 
cross section, da'^'^ [s) / dy , is defined via the eikonal (PT|) taken at the energy 




(28) 
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We use a set of parton distribution functions (PDFs) from Ref. [12] based 
on an analysis of existing deep inelastic data in the next-to-leading order 
QCD approximation in the fixed-flavor-number scheme. The extraction of 
the PDFs is performed in :16| simultaneously with the value of the strong 
coupling and high-twist contributions to structure functions. The PDFs are 
available in the region 10"^ < x < 1, 2.5 GeV^ < < 5.6 • 10^ GeV^ HE]. 
So, no extrapolation in variable x is needed. 

We put Xmin = A^/s in (j2ZI). Since the eikonal is effectively cut at 
b = bo{s) (|^. we take the mass scale in PDFs to be M = 1/boiVS). The 
effective impact parameter bo is much smaller than the size of the nucleon. 
Thus, our assumption that the neutrino interacts with the constituents of 
the nucleon and probes its inner structure, is well justified. 

The differential cross section as a function of y, the energy fraction de- 
posited from the neutrino to the nucleon, is presented in Fig. IHl for the neu- 
trino energy E^, = 10 EeV and three sets of parameters of the RS model. 

In order to estimate an effective range of variable y which contributes to 
the neutrino-nucleon cross section, we have calculated a quantity 



where yo is a minimum fraction of energy lost by the neutrino (deposited to 
the nucleon). The dependence of the quantity <Jei{y > yo) on yo at different 
values of the neutrino energy E^, is shown in Fig.[7|for = 2 TeV, fi = 0.05. 
Next two figures. Fig. |Hland Fig. ^ show the dependence of crei{y > yo) on 
the parameters /i and A^^. 

Ultra-high energy cosmic neutrinos have not yet been detected (see non- 
observation of neutrino- induced events reported by the Fly's Eye [T7], the 
AGASA and the RICE ^H] collaborations). A number of experiments 
under construction will allow to measure fluxes of such neutrinos within 
the next few years. Among them are Pierre Auger Observatory, IceCube 
neutrino telescope at the South Pole, Anita radio detector for a balloon 
flights around the South Pole, as well as EUSO, SalSA and OWL proposals. 
We will consider the first possibility |20j . 

The number of horizontal hadronic air showers with the energy Egh larger 
than a threshold energy Et^, initiated by neutrino-nucleon interactions, is 



1 




(29) 



yo 
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given by 




A{yE) 9{yE - E,,) 



+ Yl '^uAE)cT',;'^{E)A{y.E,J9{yiE-E,,)\, (30) 

i=e, fi, T 

where Na = 6.022 ■ 10^^ g-\ T is a time interval (one year, in our case), and 
A{E) is a detector acceptance as a function of a shower energy (in units of 
km^ steradian water equivalent = 10^^ g). The quantity ^^^{E) in (jHUj) is a 
flux of the neutrino of type i. Both neutrino and antineutrino are assumed 
in the sums in Eq. (j30|) . The product E^^.{E) is in units of cm~^ yr~^- 
We have taken into account that the energy of the shower resulted from the 
gravitational interaction is equal to yE, and that this interaction is universal 
for all types of neutrinos. 

For the energy distribution of the neutrino in the SM processes, we have 
used the approximation da^'^{y, E)/dy '::=l a{E) 6{y — yi). The inelasticity yi 
defines a mean fraction of the neutrino energy deposited into the shower in 
a corresponding SM process. We have put ye = I for SM charged current 
interactions initiated by electronic neutrino, while for SM neutral interactions 
initiated by Ue and for z/^/z/^-e vents we have taken ye = Vfi = Vt = 0.2 |21j . 

The number of extensive quasi-horizontal showers induced by so-called 
cosmogenic neutrinos which can be detected by the array of the southern 
site of the Pierre Auger Observatory, is presented in Table ^ for several sets 
of the RS parameters. These values of parameters are chosen in such a 
way in order not to violate experimental and theoretical bounds presented 
in Fig. 0^ The cosmogenic neutrino flux is taken from Refs. [22]; assuming 
Emax = 3 ■ 10^^ eV. The acceptance of the Auger detector is taken from 
Ref. (it is not assumed that a shower axis falls certainly in the array). 
The threshold energy Eth is chosen to be 10^^ eV. 

For comparison, a SM background is presented in the last row of the 
Table This value is in agreement with the number obtained recently for 
the same neutrino flux in Ref. [23]. 

The cosmogenic neutrino flux is the most reliable one, since it relies only 
on two assumptions: (i) the observed extremely high energy cosmic rays 

•^Remember that in our notations k/Mpi = n, mi ~ 3.83 /uA^. 
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Table 1: Yearly event rates for nearly horizontal neutrino induced showers with 
(^zenith > 70° for the cosmogenic neutrino flux from Ref. \22,^ for three sets of the 
parameters. Number of events corresponds to one side of the Auger ground array. 





A^=2 TeV, /i=0.10 


A^=3 TeV, /i=0.05 


A^=3 TeV, /i=0.10 


SM+grav 


0.81 


0.66 


0.43 


SM 


0.24 



contain protons, (ii) these cosmic rays are primarily extragalactic in origin. 
Note, however, that the cosmogenic neutrino flux may be signiflcantly de- 
pleted, if a substantial fraction of the cosmic ray primaries are heavy nuclei 
rather than protons |25j . 

The cosmogenic neutrino flux not only represents a lower limit on the 
flux of ultrahigh energy neutrinos, but it also can be used to put an upper 
limit on the neutrino flux. In Ref. [211 an upper limit (called WB bound) 
on a flux of neutrinos from compact sources which are optically thin to py 
and pp interactions (such as active galactic nuclei) has been obtained. The 
number of showers which can be registered by the Auger detector for this 
case is shown in Table |H We have chosen the same threshold energy Eth = 
10^^ eV and put E^ax = 10^^ eV. 

Table 2: The same as in Tabled but for the Waxman-Bahcall neutrino flux |26j . 





A^=2 TeV, /i=0.10 


A^=3 TeV, /i=0.05 


A^=3 TeV, /i=0.10 


SM+grav 


1.03 


0.78 


0.49 


SM 


0.28 



Note, the so-called cascade upper limit on transparent neutrino sources 
(MPR bound) is 43 times higher that the WB bound. It exploits the EGRET 
data on the diffuse gamma-ray background |28j . 

The lower bound on the cosmogenic neutrino flux was also obtained under 
assumption that the observed extremely high energy cosmic rays below 10^° 
eV are protons from uniformly distributed extragalactic sources [29J. It uses 
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the fact that the proton are accumulated around the energy E^^^^ = 4-10^^ eV 
due to the GZK mechanism [20] • The lower cosmogenic neutrino spectrum 
is practically cut at c:^ 2 ■ 10^^ eV j^H]- Other recent estimates of the 
cosmogenic neutrino fluxes can be found in Refs. I24j . 

4 Conclusions and discussions 

In the present paper we have calculated the contribution from the massive 
graviton modes to the eikonal in the RS model. The results were applied to 
the neutrino-nucleon scattering at transplanckian energies. Both differential 
and total cross sections are estimated for the different sets of the parame- 
ters of the model. By using differential cross sections, we have calculated 
the number of quasi-horizontal neutrino induced air showers which can be 
detected at the Auger Observatory per year. The estimates were obtained 
for two fluxes of cosmogenic neutrinos. 

The differential cross section, da{y)/dy, where y is the fraction of the 
neutrino energy E^ deposited to the shower, can reach tens of mb at y = 0, 
depending on energy (Fig.lH)). However, the differential cross section exhibits 
a rapid fall-off in y, starting at some small y. As a result, the gravitational 
cross section appears to be approximately one order of magnitude larger than 
the SM cross section at the same energy. To illustrate this statement, let us 
fix the parameters of the RS model to be A^ = 2 TeV, fi = 0.1. Then we 
have {da{y)/dy)\y=o ~ 4 mb for E^ = 10^° GeV (see Fig. E}. As one can 
see in Fig. El da{y)/dy begins to fall rapidly at y > 10^^. The numerical 
calculations show that a ~ 4 ■ 10~^ mb for this case (dashed curve in Fig. |H)). 

The energy of neutrino induced air shower, Egh = yE^, is bounded from 
below by a threshold energy Et^. Thus, the fraction y should obey the in- 
equality y ^ Eth/Emax, where Emax is a maximum energy in a neutrino 
spectrum. For Eth = 10^ GeV and Emax = lO^^^^^) q^v, we get y ^ IQ-^^^l 
Thus, the air showers event rate is defined by the region of y, in which 
neutrino-nucleon cross section da{y)/dy is significantly reduced in compari- 
son with its magnitude at y = 0. Nevertheless, gravity contribution to the 
event rate at the Auger detector is several times larger than the SM back- 
ground, as one can see from Tables [H and El 

Recently, model independent bounds on the inelastic neutrino-nucleon 
cross section derived from the AGASA ^H] and RICE search results 
on neutrino events were obtained J24j. The bounds exploit the cosmogenic 
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neutrino fluxes from Refs. [221 121] ■ However, tliey were derived under an 
assumption tliat tlie total neutrino energy goes into sliower energy, tfiat is 
y = 1. As we liave seen, it is not a case for the gravitational interactions 
originated from t-channel KK gravitons, which prefer y -C l."^ Generally, in 
order to extract an upper limit on atot, a dependence of da{y)/dy on y is 
needed. So, we conclude that the bounds from Ref. [21] can not be directly 
apply to the neutrino-nucleon cross sections derived in our scheme. 
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Appendix 

In Appendix we calculate the dependence of the eikonal on variable 
u firUj) . Let us define 

oo 

I{u) = j dxJo{ux) R{x), (A.l) 



with 

Rix) = ^ (A.2) 

being the ratio of two modified Bessel functions. 

It easily to see from ()A.1|) that I{u) —* u^^ at m ^ 0, since R{x) — » 1 at 
X — s> oo. The asymptotics of I{u) at large u is defined by the behavior of the 
integrand at small x which looks like 

^ 7 

+ ^ ^ + 0(x9). (A.3) 

^ ^ 4 96 1536 23040 ^ ' ^ ' 

Let us now demonstrate that at m — oo the function l{u) (jA.ll) decreases 
faster that any fixed power of m^^. By using well-known relation j22j 

x^-V,_i(Ma;) = - \x^JAux)\ , (A.4) 
u \xdx/ 

^In processes initiated by graviton t-channel exchanges in large extra dimensions mean 

energy loss is also small, as was pointed out in Ref. [HI]. On the contrary, in a process of 

black hole production, the neutrino loses most of its initial energy (y ~ 1). 
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and integrating ()A.1|) by parts k times, we obtain: 

oo 

I{u) = ^ j dx Jk{ux)Fk{: 



where Jk{z) is the Bessel function, and 

Fk{x) = {-If x^'^' 



(A.5) 



d 



R{x) 



X 



(A.6) 



The function Fk{x) ()A.6|) has the following properties: it is proportional 
to x'''^^ at X — i> 0, and it decreases as x~'' at x — oo. For any positive 
integer k, it depends only on x and on the ratio R{x) ()A.2|) (but not on Ii{x) 
and hix) separately) due to the following relations between modified Bessel 
functions 



d_ 

dx 



hix) 



l2{x) + -h{x) 

X 

h{x)--h{x) 

X 



For instance, for k = 1 one has 
Fiix) = 



X 



(A.7) 



(A.8) 



This expression has asymptotics a;^/48 and x ^ at small and large x, respec- 
tively. For k = 2 one gets 

4i?(x) 



F.ix) 



6 



1 - 



X 



- 2x R{x) + 12 [R{x)] ^ + 2x [R{x)] ^ | (A.9) 



The asymptotics of F2{x) are x^/192 and 3x~^. 

Since k is an arbitrary positive integer, we conclude from ()A.5|1 and ()A.6|1 
that lim I{u) = for any a > 0. 

u—*oo 

The integral in ()A.5|1 , contrary to an original one (jA.ljl , converges rapidly 
at X — >^ oo for k ^ 2, and could be used for numerical calculations. It 
cannot be calculated analytically. However, there exists an expression which 
approximates our integral with a very high accuracy: 



I{u) = J dx Jo{ux) R{x) , 





(A.IO) 
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with 

The function R{x) has the following expansion at < 12 (compare with 
Eq. da): 

^ ^ _ ^ + ^ ^ + 0(x9). (A.12) 

^ ^ 4 96 1536 110592 ^ ' ^ ' 

The integral ()A.10|) is a table one 133]: 

/(m) = — exp(-2y3M). (A.13) 
2u 

We can integrate the RHS of ()A.10|) twice by parts, 

oo 

I{u) = ^ f dxh{ux)F2{x), (A. 14) 



and compare F2{x) = (3-\/3/2) x^ jix?' + 12)^/^ with the corresponding func- 
tion F2{x) ()A.9j) . The result of our calculations is presented in Fig. Uni 
Formula ()A.13|) gives practically the same dependence on variable m as a 
numerical integration of the exact expression by using formula ()A.5|) (with 
k = 2) does, see Fig. [TT]^ Thus, I{u) exhibits an exponential fall-off (as 
we expected, see above), and it becomes as small as I{u) ~ 0.01 already at 
u = 1.2. 

Taking all said above into account, we put I{u) — > I{u), that results in 
the analytical expression for the eikonal presented in the text (PT|) . 
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Figure 1: Experimental and theoretical constrains on the RS model in the 
two-parameter plane n/Mpi - mi [3,. The allowed region lies in the center 
as indicated. 
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s = 2 10" GeV 
s = 2 1 0'° GeV 
s = 2 1 0' GeV' 




Figure 2: The differential cross section as a function of dimensionless variable 
y for three fixed values of the invariant energy. The parameters of the RS 
model are chosen to be = 2 TeV, and /i = 0.05. 



A, = I TeV 
A, = 2 TeV 
A = 4 TeV 




Figure 3: The differential cross section as a function of variable y for three 
values of the mass scale at fixed energy s = 2-10^^ GeV^ (with the 
parameter /x = 0.1). 



19 




Figure 4: The differential cross section as a function of variable y for three 
values of the RS parameter /i at fixed energy s = 2 ■ 10^° GeV^ (with the 
scale = 1 TeV). 




Figure 5: The reduced (dimensionless) differential cross section as a function 
of variable y for different parameter sets (A^, /x) at invariant energy s = 2-10^° 
GeV^. The product /lA^ is taken to be 100 GeV (200 GeV) for three first 
(two last) sets. 
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Figure 6: The differential neutrino-proton cross section as a function of 
tfie fraction of the neutrino energy deposited to the proton. 
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Figure 7: The inelastic neutrino-proton cross-section obtained by integrating 
the differential cross section in the region ?/o ^ 2/ ^ 1 as a function of yo, the 
minimal fraction of the neutrino energy deposited to the proton. 
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Figure 8: The same as in Fig. [7|but for the different value of yU = 0.1. 



-2 



10 



10 







100 EeV 


3 


-- c = 


10 EeV 


.... E, = 


1 EeV 


4 












5 


! A„ = 4 TeV. = 0.05 




6 








-6 -4 -3 -2 

10 10 10 10 10 10 



Yo 

Figure 9: The same as in Fig. [7| but for the different value of = 4 TeV. 
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Figure 10: The exact integrand vs. approximate one (after the integration 
of both integrals by parts twice). See Appendix for details. 




Figure 11: The input integral vs. approximate one as a function of the 
dimensionless variable u (fT^ . See Appendix for details. 
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